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ABSTRACT 

Context. Variations in the observed dust emission and extinction indicate a systematic evolution of grain properties in the transition 
from the diffuse interstellar medium (ISM) to denser molecular clouds. 

Aims. The differences in the dust spectral energy distribution (SED) observed from the diffuse ISM to denser regions, namely an 
increase in the spectral index at long wavelengths, an increase in the FIR opacity, and a decrease in temperature, are usually assumed 
to be the result of changes in dust properties. We investigate if evolutionary processes, such as coagulation and accretion, are able to 
change the dust properties of grains in a way that is consistent with observations. 

Methods. We use a core-mantle grain model to describe diffuse ISM-type grains, and using DDA we calculate how the accretion of 
mantles and coagulation into aggregates vary the grain optical properties. We calculate the dust SED and extinction using DustEM 
and the radiative transfer code CRT. 

Results. We show that the accretion of an aliphatic carbon mantle on diffuse ISM-type dust leads to an increase in the FIR opacity by 
a factor of about 2 and in the FIR/submm spectral index from 1.5 to 1.8, and to a decrease in the temperature by about 2 K. We also 
show that the coagulation of these grains into aggregates further decreases the temperature by 3 K and increases the spectral index up 
to a value of ~2. The FIR opacity is increased by a factor of 3 (7) for these aggregates (with an additional ice-mantle) compared to 
the diffuse ISM-dust. 

Conclusions. Dust evolution in the ISM resulting from coagulation and accretion, leads to significant changes in the optical properties 
of the grains that can explain the observed variations in the dust SED in the transition from the diffuse ISM to denser regions. 

Key words. ISM: dust, extinction, ISM: evolution, ISM: abundances 


1. Introduction 

Indications of dust evolution from the diffuse interstellar 
medium (ISM) to molecular cloud environments (n > 10^ 
H/cm^) have rece ntly been ob served with the H erschel 
dF'ilbratt et al.l l201()ll and Planck (iLamarre et al.l l2010l) space 
observatories in the far infrared (FIR) and submillimetre 
(submm) wavelength range. Earlier, dust evolution was al¬ 
ready indicated with ISO, Spitzer, SPM/PRONAOS, and IRAS 
data. For example, a temperature decrease of large grains, 
which are in thermal equilibrium with the radiation held, 
from 2 0.3 K in the diffuse ISM dPlanck Colla boration et al.l 
2014al) to less than 14 K in denser regions (e.g., Lagacheet^ 


1998 : IStennik et al.l 12003^ iPlanck Collaboration et al.l201 Ibh 


was observed together with an increase in the dust opacity 
by a f actor of 2 to 4 in the FIR/submm (e.g., Stepnik et al.l 
2003|: I Planck Collaboration et al.l l201 Ibt Lluvelaet^l 1201 


Rov et al.n2013 ). Additionally, the mid-IR emission, which 


vironment (Dupac et al. 20031; Stepnik et al 

I 2 OO 3 I: Desert et al. 

2008; 

iParadis et al.l 1201 Ol Veneziani et al.l 

2 OIOI; Bracco et al. 

2011; 

Planck Collaboration et alJ 12011 hi 

JuvelaetalJ 120131; 

Sadavov et al.ll2013l). Values for B and T are derived bv htting 


come s mainly from stochastically heated very s mall grains 
(e.g. iDraine &Xil 1200 ll ICompi^gne et al.l 1201 ll: Llones et ^ 
I 2 OI 3 I) . de creases by 80-100% from diffuse to denser re¬ 
gions (e.g.. lLaureiis et al.lll991l : [^pnik et alJ2003l:ljuvela et al.l 
I 2 OI 2 I) . In many of these papers, an anti-correlation of the 
dust colour temperature T with the opacity spectral index [5 


observations at wavelengths longer tha n around 60 um wi th a 
modihed blackbody. As was shown by lJuvela et~an (120131) the 
values derived for j3 and T depend strongly on the htting method 
and the wavelength range of the observations where the modi¬ 
hed blackbody has been htted. In general, the;^^^ approach over¬ 
estimates p and underestimates T when htting noisy observa¬ 
tions. In the lite rature, we therefor e hnd a wide range for value s 
of p, from 0.8 (iDuirac et al. 20()3l) to 5 (IVenez iani et al j2010l) . 
an d of r, from 7^ (iDesert et al.ll200^ iParadis et al.ll201ol) to 80 
K (iDupac et al.ll2003l) . A systematic study of the various htting 
methods (x^, hierarchical models , and the Bayesian method) was 
carried out bv Lluvela et al.l (l2013l) . Their main result was that the 
observed anti-correlation cannot be explained only by bias in 
the htting methods but must originate, at least partly, in intrinsic 
variations in dust properties. It is usually assumed that all four 
changes in the observed SEDs occur because of the evolution of 
dust from diffuse to dense regions. These variations result from 
the evolutionary processes acting on dust as it responds to its 
environment. 

Changes in the environment include variations in the lo¬ 
cal gas density and in the radiation held. Coagulation, accre¬ 
tion, fragmentation, sputtering, and photoevaporation change 
the properties of dust grains, such as their structure, shape, 
material composition, and size. This results in a modihca- 
tion of the optical properties of the grains and therefore to 
changes in the spectral energy distribution (SED). Variations 
in dust optical properties resulting f rom coagulation have 
been explored i n many papers (e.g. iBazell & Dw^ ]_990t 
IOssenk(^ll993HOssenkopf & Henninglll99^IStognienko et alJ 
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1995 


2011 


Dwek 1997 


I Kohler et al 


, Fogel & Leund Il998l : lOrmel et aP l2009l 
201 1[ 2012i) . In the transition from the dif¬ 


fuse ISM to denser regions, the main evolutionary processes 
are the accretion of mantles composed mainly of amorphous 
carbon, water ice, and other molecular species (Jones etalJ 


19^ lMathis|[l9^ IWhittet et al.lll996l: Idel Burgo ^ 
Kiss_eLair2QQ 0 and grain coagulat i on jB epiardeLal 


del Burgo et alJ 1255^ Stepnik et alJ 20031: iKiss etTlI 


2003 


1999; 


2006t 


Ridderstad et alJ 200d iParadis et alJ 120091 1 Ysard et alj|2013l) . 


Observational evidence o f grain coagulation in denser clouds 
was found by lJural (Il980ll through an analysis of the extinction 
towards p Oph. 

The yS-T anti-correlation is seen for amorpho us carbon and 


silica t e materials in laboratory measurement s (|Agladze_eiM 


1996 ; iMennella et al.ll 19^ iBoudet et alJ 2005t ICoupeaud et~ 



201 ih and the DCD/TLS model dMenv et al. l2007h . We do not 
take variations of j3 with grain temperatures into account in our 
model, but plan to in future studies. Since these variations occur 
mainly at long wavelengths and for large temperature variations, 
the influence on the SED and extinction are likely to be small for 
our assumed dust composition and size distribution. 

However, so far, no dust model or dust-evolution model have 
been able to self-consistently explain all four observed changes 
in the SED from the diffuse ISM to denser regions, i.e. changes 
in temperature, spectral index (and their anti-correlation), opac¬ 
ity (emissivity) and the mid-infrared emission. Model calcula¬ 
tions taking into account the coagulation of separate grains into 
aggregates could explain the decrease in temperature and the in¬ 
crease in opacit y, but were not able to explain the observed in¬ 
crease in p (e.g. iKohler et alj|2012|) . This theoretical study was 
later confi rmed by simulatio ns including radiative transfer cal¬ 
culations (lYsard et alj|2012h and the analysis of Herschel ob¬ 
servations of_a_dei^ filument in the Taurus molecular cloud 
complex (lYsardetalJ 2013h. T hese three studies were based 
on the ICompiegne et al. (1201 Ih dust model, which consists of 
neutral and charged PAHs, very small grains (VSGs, ~ 1. 2-15 
nm in radius) of amorphous carbon (IZubko et al.l Il996h and 
big grains (BGs, 15- 110 nm in radius) of astrono mical silicate 
( Draine & Lee! Il984h and of amorphous carbon (IZubko et al.l 
119961) . 

Recently, new opt ical constants for hydrogenated amorphous 
carbon were derived (Ijonesll2012allbll3l . A new dust model that 
uses these materials, in conjunction with new silicate data, ex¬ 
plains the observational data of dust in the di ffuse ISM, such as 
the IR to EUV extinction and the dust SED (I.Iones et al.ll2013t 
iKohler et al.l 120141 : lYsard et al.ll2015h . In the present study, we 
use this new dust model and coagu late the grains into a ggregates 
following the same approach as in lKdhler et al.l (l2012h . We also 
consider the accretion of carbonaceous material and ice from the 
gas phase onto the surfaces of the aggregates. 

This paper is organised as follows: In Sec.|2]we present the 
dust model and dust evolution. In Sec.[3]we describe the method 
and calculations and present the results in Sec.|4] We summarise 
and conclude our study in Sec.|5] 


2. Dust model and dust evolution 

Our study is base d on the recently published dust model from 
iJones et d] (120131) for the diffuse ISM («h < 100 H/cm^). In 
the diffuse ISM, we assume that VSGs and BGs are mainly sep¬ 
arated. Very small grains up to about 20 nm in radius consist 
purely of aromatic-rich amorphous carbon because of UV photo¬ 
processing. Larger VSGs and BGs of amorphous carbon show a 


Size in nm 

Fi g, 1. The mass dist ribution of the CM dust grains as described 
in ijones et al.l (1201 3h . of core-mantle-mantle (CMM) grains and 
of aggregates (AMM) and aggregates with an accreted ice man¬ 
tle (AMMI). 


core-mantle structure, where the mantle with a thickness of ~20 
nm is aromatic rich, while the core is aliphatic rich. Another BG 
population consists of amorphous silicates, with Ee/EeS nano¬ 
particle inclusions and a mantle of aromatic-rich amorphous car¬ 
bon, which is formed by accretion and the coagulation of VSGs 
and has a thickness of ~5 nm. We assume amorphous silicates 
of forsterite-type and of enstatite-type and the inclusions occupy 
a volume of 10% and consist of a mix ture of 30% EeS and 70% 
Ee as described in IKohler et alJ(l2014l) . Starting from these grain 
populations, we consider dust evolution based on the processes 
of accretion and coagulation, which we assume to occur in rela- 
tively dense inters tellar environments (hh > 1500 H/cm^, as in 
lYsard etalJdloIl l. 


2.1. Accretion of carbonaceous mantles 

The first process that we consider is accretion. Because of the 
accretion of C and H atoms from the gas phase, an additional 
amorphous carbon mantle is formed on the surface of grains 
(I.Iones et al.ll2014l) . IParvathi et al.l (l2012l) sho w that more car¬ 
bon i s available than usually assumed (see e.g. lComniegne et’aP 
IMTl) . so that there appears to be sufficient gas phase carbon to 
form such a carbonaceous mantle. In low-density regions, this 
newly accreted m antle would arom atise quickly because of UV 
photo-processing (iJones et al.ll2014l) . As a result, the grain prop¬ 
erties would not change dramatically, since the aromatic amor¬ 
phous mantle would only increase the particle size but would not 
change the material properties. This would have the strongest in¬ 
fluence for the small grains, where the mantle would increase the 
particle si ze by a larger fact or. Howeyer, this case was already 
studied in lYsard et al.l(l2015l) . In denser regions, where UV pho¬ 
tons are significantly attenuated (Ay > 1), the accretion of C 
and H atoms will result i n an aliphatic-rich amorphous carbon 
mantle (IJones et alJl2014t) . so that BGs will haye a core-mantle- 
mantle structure. We inyestigate how the optical properties of the 
grains change because of the accretion of such an aliphatic-rich 
amorphous carbon mantle. 

In the following, we use abbreyiations to describe the parti¬ 
cles with accreted mantles: 

CM: Core-mantle grains from the IJones et ^ (1201 3l) and 
IKohler et al.l (1201 4l) model. 
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Fig. 2. 2abs (solid lines) and Qsca (dashed lines) for CMM grains with cores of amorphous carbon (left) and cores of amorphous 
silicate of olivine type (middle) and pyroxene type (right) for different grain radii in yttm. The inner mantle is aromatic-rich and the 
outer mantle aliphatic-rich amorphous carbon. At long wavelengths Qabs and Qsca increase with increasing grain size. 



Fig. 3. 2abs (solid lines) and Qsca (dashed lines) of AMM structure (left) and of AMMI structure (right) for different radii in yum of 
volume-equivalent spheres ay At long wavelengths Qabs and Qsca increase with increasing grain size. 


CMM: Core-mantle-mantle grains, which are the original grains 
with an additional aliphatic-rich mantle. 


2.2. Coagulation into aggregates 

In denser regions of the ISM (mh > 10^ H/cm^, Ay > 1), 
grains are expected to coagulate into aggr egates. Considering 
the time scales (see e.g. iKohler et al.ll201^ . the coagulation of 
VSGs onto the surface of BGs is relatively fast, so that the man¬ 
tles on the BGs thicken with time. We assume that all VSGs 
coagulate onto the surfaces of BGs and that the population of 
VSGs is therefore incorporated into the BG population aggre¬ 
gates. The time scales t o coagulate a few BGs into aggregates 
(f = 4 X 10^ - 10^ yr, iKohler et al.ll2^ are consistent with 
typical cloud lifetimes (f ss 10^ years. IWalmslev et al.ll2004ll . 
Furthermore, amorphous carbon mantles tend to favour the co¬ 
agula tion of BGs, since they act as a sort of glue (IChokshi et al.l 
Il993h . Observations indicate that the relative velocities between 
the BGs can be small and hence fragment ation processes mini¬ 
mal (lOrmel et al.ll2009HPagani et akOOlOh . 

We assume aggregates that are formed of 4 BGs (monomers) 
in an i ntermediate elongated shape, as described by IKohler et al.l 
(l2Q12h . At long wavelengths, the difference in optical properties 
compared to separated grains, is largest for the coagulation of a 
few mon omers and remain s rather constant for larger aggregates 
(see e.g. lMackowskill^OOht IKohler et al.ll201lh . The time scales 


to coagulate 4 BGs are realistic in comparison to the cloud life- 
ti me in relatively dens e media {tin - 4 x 10^ cm“^) as shown 
by IKohler et al.l (l2012h . With the increasing density of the en¬ 
vironment, the coagulation rate between BGs increases so that 
larger aggregates are able to form. Larger aggregates have lower 
temperatures, but the differences in /3 and th e FIR opacity are 
small compared to smaller aggregates (see e.g. lMackowskil20()ra 
IKohler et al.l2011|) . We therefore only consider small aggregates 
consisting of 4 BGs in this study. 


The constituents of the aggregates are one BG of aliphatic- 
rich amorphous carbon and three BGs of amorphous silicate of 
olivine and pyroxene type with Fe/FeS inclusions as described 
above. This distribution a grees with the ma ss fractions of BGs 
in the ISM as derived by Jones et ^ (l20I3l) . We then consider 
the accretion of a second mantle in denser regions as in Sec. 13.II 
This second mantle consists of aliphatic-rich amorphous carbon, 
which is unprocessed because of the high UV extinction. 


W ith increasing density (line of sight Ay > 1.5. I.Iones et al.l 
l2014t) . the formation of a water-ice mantle occurs. We therefore 
assume the formation of such ice mantles on the surface of the 
grains after forming the aggregate. 


In the following we use abbreviations to describe the coagu¬ 
lated aggregates: 

AMM: Aggregates consisting of CMM grains. 

AMMI: Aggregates consisting of CMM grains with an addi¬ 
tional ice mantle. 
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Fig. 4. The differences in and Qsculctv for the evolution of grains with CM grain radii of 0.06 fjm. Left: (dashed- 

dotted and solid lines) and Q^^^lay (dotted and dashed lines) for CM and CMM grains, respectively. Right: (2abs/av (solid and 
dashed-dotted lines) and Qsca/ay (dashed and dotted lines) for CMM grains (as in left pannel) and for AMM and AMMI grains, 
respectively. 01 indicates amorphous silicate of olivine-type, Py indicates amorphous silicate of pyroxene-type and Ca indicates 
amorphous carbon. 


3. Theory and calculations 

Model calculations to derive the optical properties of aggregates 
are carried out with the disc r ete-dipo l e app r oximation (PDA 
Purcell & Pennvnackerl 119731 : iDrainel Il98^ iDraine & Flataul 


20101) . We u se the PDA code DDSCAT 7.1 0 (for a detailed de¬ 
scription see lDrain^ll988tlbraine & Flataul2 010() including the 
lattice dispersion relations (LDR) ( Draine & Goodmanll993[) . 


3.1. Accretion 

In order to study the accretion process, we consider grain radii 
from 0.04 to 1.0 /urn. Up to a radius of 0.441 pm we use DDA 
to derive the optical properties. For DDA calculations, the grains 
are represented by interacting dipoles. In order to obtain the cor¬ 
rect size ratios between BGs and VSGs, we assume that BGs 
of radius (cores) 0.02 pm are 7 dipoles across, BGs of 0.06 pm 
are 17 dipoles across and all larger BGs are 35 dipoles across. 
Since we want to coagulate these grains into aggregates contain¬ 
ing 4 BGs, we cannot consider more dipoles because the DDA 
calculations are then computationally too expensive. We assume 
that 1000 dipoles, 3 nm in radius, form the original mantle plus 
the first coagulated VSGs of aromatic amorphous carbon and 
that another 1000 dipoles of aliphatic amorphous carbon, 3 nm 
in radius, form the accreted mantle on each BG. For cores of 
0.02, 0.06 and 0.11 pm, where each dipole is 3 nm in size, 1000 
dipoles are distributed on the grain surface. For larger BGs, the 
dipole size is larger than 3 nm, since we keep the number of 
dipoles across the BG constant but increase the radius. We there¬ 
fore reduce the number of dipoles in the mantle when increasing 
the grain size, so that the mantle mass is constant for each BG . In 
order to avoid an antenna effect (see e.g. lKohler et al.ll2012l) . the 
mantles form a layer so that chains of dipoles protruding from 
the surface are prevented. The optical property re sults are av¬ 
eraged over 125 grain orientations as discussed in iKohler et al.l 
(l2Q12l) . which allows a relatively fast calculation and gives suf¬ 
ficient accuracy. 

For grains with radii larger than 0.4 41 pm, we use the 
Maxwell-Garnett effective medium theory (iBohren & Huffm^ 


Il983l) to derive the best mixture of aromatic and aliphatic amor¬ 
phous carbon in the m antle and amorphous silica te in the core. 
We then use BHMIE (iBohren & Huffm^ll983l) to derive the 
optical properties for the grainfl The best mixture was derived 
by comparing the optical properties for smaller grains, where 
DDA results are available. We find the best mantle mixture of 
50% aromatic-rich and 50% aliphatic-rich amorphous carbon in 
volume and a mantle to core volume relation of 0.08% to 99.2%. 

The mass distribution of CMM grains is shown in Fig.[T] The 
range of the mass distribution is determined by a^nm, which is the 
smallest grain size which is possible when adding the accreted 
mantle and by t/max, which is the largest grain size we calculate. 
We see clearly an increase in the mass as a result of the accre¬ 
tion of material. For the CMM grains, we use 182 ppm of C for 
the first mantle and 224 ppm of C for the second mantle. The 
total o f 406 ppm is in good agreement with the c arbon abun¬ 
dance (ICompi^gne et al.ll201 UiParvathi et ^l2012l) . We derived 
the average density to be 1.46 g/cm^ for a core of amorphous 
carbon and and 1.59 g/cm^ for a core of amorphous silicate. 


3.2. Coagulation 

We now coagulate the CMM grains as described in Sec. l3.1l into 
aggregates of 4 BGs, where the cores of 2 BGs are of amorphous 
silicate of olivine type, 1 BG of amorphous silicate of pyrox¬ 
ene type and 1 BG of amorphous carbon. The optical properties 
of the aggregates are calculated with DDA for six sizes, with 
volume-equivalent sphere radii, ay, of 0.05, 0.06, 0.09, 0.17, 
0.318, 0.476 and 0.70 gm, averaged over 125 aggregate orien¬ 
tations. The aggregate size distribution is shown in Fig. [1] The 
maximum is shifted to larger grains because of the coagulation 
process. With our adopted formalism the smallest aggregates, 
with radius flmin, actually only consist of carbonaceous mantle 
materials. We assume amax = 0.70 pm, since it was computa- 

* Using BHCOAT dBohren & Huffmaiil 1 198.31) is not possible, since 
the mantle on the large grains would be too thin to derive reasonable 
results. 
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Fig. 5. The SED of the CM grains, CMM grains and aggregates, 
AMM and AMMI. 


tionally too expensive to calculate the optical properties of larger 
aggregates. We derive the averaged density to be 1.55 g/cm^. 

In a second step, we also study the changes in the optical 
properties when accreting an ice mantle on the surface of the 
aggregates. For the ice mantles we assume 6000 dipoles of 3 nm 
in size per BG. We take the optical constants of ice from lWarrenI 
(Il984l) . The mass distribution for these icy aggregates is also 
shown in Fig. [1] The a^m is again shifted to larger grain sizes. 
The increase in mass due to the accreted ice mantle is clearly 
visible. The average density decreases to 1.18 g/cm^. For the ice 
mantle we use 500 ppm of O which is in a greement with the 
abundances give n inIComp iegne et al.l (1201 Ih . 

As shown bv iKohler et al.l (1201 ih . the optical properties of an 
aggregate depend on the connections between the monomers. In 
this study we consider a contact area of 3, 5 and 9 dipoles across 
forBGs of 7,17 and 35 dipoles across, respectively. We also con¬ 
sider that the BGs form the aggregate and VSGs are connected 
to their surf aces, so that in the contact area the BGs “touch”. 
As shown bv lKdhler et al.l (l201 ih . the increase in the absorption 
coefficient depends on the real part, n, of the optical constants. 
This is at the origin of the increase in opacity (emissivity) at 
250 fim. Since the amorphous silicate and the amorphous car¬ 
bon (Eg=0.1) have similar values of n at 250 /rm, the differences 
in the absorption coefficient are small when VSGs provide the 
contact area. This was tested by using dimers. Initially, we con¬ 
nect silicate BGs and coagulate VSGs of amorphous carbon on 
their surfaces. In follow-up calculations, we separate the silicate 
BGs by two dipoles and fill the contact area and the surface of 
the BGs with VSGs of amorphous carbon. The number of VSGs 
is the same in each case. The absorption coefficients vary at most 
by only 4% and 5% for aromatic and aliphatic amorphous car¬ 
bon, respectively. 


3.3. Dust emission 

Firstly, we use the DustEM cod^^ (IComniegne et al.l 1201 Ih to 
calculate the SEDs and the extinction of our grains. We assume 
a diffuse (transparent) medium illuminated by the interstellar ra¬ 
diation field (ISRF) of iMathis et al.l (Il983h for a galacto-centric 
distance of 10 kpc. 

^ DustEM is a publicly available numerical tool to calculate the SED, 
the extinction and the grain temperature distribution for dust in the ISM, 
available at http://www.ias.u-psud.fr/DUSTEM/ 




Fitting intervai 

Fig. 6. Colour tempe rature, spectral ind ex 0 and opacity derived 
with DustEM for the lJones et akl (12013h and iKohler et alJ (l2014l) 
model (CM), CMM grains and AMM and AMMI aggregates. 


Secondly, we simulate denser environments by coupling 
DustEM to the 3D contin uu m radiative transfer code CRT as 
described in iJuvelal (l2005h and lYsard et al.l (l2012h . We consider 
spherical clouds with the following radial hydrogen density dis¬ 
tribution, 

”»('■) = ■;— ,2 ^ 

1 + (t'/Rflat)2 

where nc is the central density, Rfiat is the internal flat radius, and 
Rout the external radius set to 0.5 pc. We define the (nc^^flat)- 
parameters for the clouds to be at equilibrium, with thei r masses 
per unit length equal to the critical mass defined by lOstrikerl 


5 
































































Kohler et al.: Dust evolution in the transition towards the denser ISM: impact on dust temperature, opacity, and spectral index 


Table 1. The temperature, spectral index, yS, and increase in opacity derived with the modified blackbody approach, e®®, at 250 /tm 
for the separate grains and aggregates, derived in different wavelength ranges. 



[100,500] 

[100,870] 

Herschel PACS & 
IRAS bands 

Planck HFI & 
SPIRE bands 


T[K] 

P 

^BB 

Mnc 

T[K] 

P 

^BB 

^inc 

T[K] 

P 

^BB 

Mnc 

T[K] 

R 6 ®® 

me 

CM 

20.2 

1.5 

1 

20.2 

1.5 

1 

20.3 

1.5 

1 

20.9 

1.4 1 

CMM 

18.5 

1.8 

2.3 

18.5 

1.9 

2.3 

18.8 

1.8 

2.3 

18.9 

1.8 2.4 

AMM 

17.5 

1.9 

2.9 

17.4 

1.9 

2.9 

17.5 

1.9 

2.9 

17.5 

1.9 3.2 

AMMI 

16.7 

2.0 

7.0 

17.3 

1.8 

7.0 

16.7 

2.0 

7.0 

17.2 

1.8 7.0 


(Il964ll . Merit = 2 C 5 /G, where G is the gravitational constant. The 
sound speed cs is computed for a gas temperature of 12 K. These 
spherical clouds are then embedded in a more diffuse medium 
with 100 H/cm^. We simulate a series of clouds with central vi¬ 
sual extinctions of Ay ~ 0.1 to 20 and illuminate them with the 
standard ISRF. 

In order to study the variations in dust populations inside the 
clouds, we consider the following models: 

Model CM: No variation in the dust properties and the diffuse 
medium as well as the s pherical clouds are filled with CM d ust 
of the diffuse-ISM type J.Iones et^l2013HKohler et al.ll2014t) . 
Model CMM: The diffuse medium is filled with CM dust while 
the spherical clouds are filled with CMM dust. 

Model AMM: The diffuse medium is filled with CM dust, and 
the dense cloud with CMM dust for Uh <1500 H/cm^ and AMM 
for higher densities. 

Model AMMI: Same as Model AMM, with the central region 
filled with AMMI grains. 

To compare our dust models with observational results, we 
finally fit our modelled SEDs with a modified blackbody, 

Sy CC V^By(T), (2) 

where By is the Planck function, p the spectral index, and T the 
colour temperature. We perform the fits in the Planck-HFl bands 
at 857, 545, 353, and 143 GHz, and the IRAS 100 jum filter. 
We also make fits to the Herschel PACS and SPIRE filters at 
100, 160, 250, 350, and 500 jim. A last pair of fits is performed, 
without instrumental filters, continuously between 100 and 500 
/rm, or 100 and 870 jim. We do standard weighted least square 
fits, where all bands/wavelengths are given an equal weight (i.e. 
we do not consider the effect of noise), in order to determine T 
andyS. 

4. Results of model calculations 

4 . 1 . Absorption and scattering coefficients 

In Eig. |2]the calculated absorption coefficients Qabs and scatter¬ 
ing coefficients Qsca for the CMM grains with different sizes and 
different core materials (amorphous carbon left, amorphous sili¬ 
cate olivine type middle and pyroxene type right) are presented. 
A difference in and Qsca with grain size is clearly visible for 
all core materials. The peak in Qsca shifts to longer wavelengths 
with increasing particle size. Qabs increases with increasing grain 
size for amorphous silicate cores and decreases at wavelengths 
longer than 1 yum for the amorphous carbon core. This decrease 
is due to the more dominant core material when the grain sizes 
increase. Eor amorphous carbon cores an increase in the fea¬ 
tures at 3.4 and around 7 yum can be seen. Eor grains with an 
amorphous silicate core, an increase in the strength of the sili¬ 


cate features can be observed and differences in the slope at long 
wavelengths are evident (A > 50 yum). 

In Eig. [3] 2abs and Qsca are presented for different sized ag¬ 
gregates, AMM (left) and AMMI (right). With increasing grain 
size both Qabs and Qsca increase. With increasing grain size the 
silicate features become more pronounced. The differences in 
the slope at long wavelengths are strongest for small aggre¬ 
gates. The smallest aggregates only consists of amorphous car¬ 
bon material and therefore show strong hydrocarbon features in 
the spectrum. Eor AMMI grains, the ice features are clearly visi¬ 
ble for small aggregates and become less evident with increasing 
grain size while the silicate features become more pronounced. 
The smallest grains consist of the pure carbonaceous AMM with 
a thick ice mantle and we find that the ice features clearly domi¬ 
nate the spectrum. 

In Eig. |4] we present the variation in the optical properties 
due to the evolution of CM grains with a radius of a = 0.06 yum. 
It can be seen that Qabs and Qsca vary with accretion (left) and 
coagulation (right). The accretion of an aliphatic amorphous car¬ 
bon mantle on the amorphous carbon grains shows a decrease in 
( 2 abs and a change in slope at long wavelengths, while differences 
in 2sca are rather small. Eor grains with an amorphous silicate 
core, we find a decrease in the silicate feature and an increase 
between 5 and 10 pm in Qabs due to the accretion of carbona¬ 
ceous material. Eor AMM grains, Qabs and Qsca increase com¬ 
pared to the single grains, but the overall shape of the curve is 
comparable to the curves of the amorphous silicate CMM grains. 
Eor AMMI grains, the ice features are clearly visible in Qabs- 
Compared to AMM grains, Qabs decreases between 0.1 and ~5 
yum and for wavelength larger than 200 pm. Eor Qsca the differ¬ 
ences are small compared to the AMM grains. 

4.2. Temperature, spectrai index, and opacity derived with 
modified biackbody fits 

Using the optical properties described in the previous sec¬ 
tion, we calculate the corre sponding SEDs fo r grains illumi- 
nated by the stand ard ISRF (iMathis et al.lll983h using DustEM 
(IComniegne et aHl201 ill . The results are presented in Eig. |5] for 
CM, CMM, AMM and AMMI grains. The SED of the CM grains 
is the one presented in iKohler et al.l (l2014l) . where VSGs are still 
separated and contribute to the emission up to ~70 pm. At longer 
wavelengths, the SED of CMM grains shows an increase in the 
intensity of the far-IR emission and a shift to shorter wavelengths 
indicating an increase in temperature due to stronger absorption 
at short wavelengths compared to CM grains. The steepening of 
the slope (i.e. an increase in the EIR opacity spectral index) is 
also clearly visible. 

Eor AMM grains the SED maximum is shifted to longer 
wavelengths, which is an indication of the decrease in temper- 
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Fig. 7. Infrared extinction of single grains and aggregates normalised at 10 fim (left) and at the K-band (right) compared to obser¬ 
vational data. The normalisation wavelength is indicated by the star. Solid lines: CM grains (black), CMM grains (green), AMM 
grains (red), AMMI grains (blue). 


ature when coagulating grains. We can also see a small bump at 
shorter wavelengths (magenta curve) which is due to the smallest 
grains which only consist of carbonaceous material. The temper¬ 
ature of these grains is around 52 K while all other AMM aggre¬ 
gates are colder than 20 K. Because of their high temperatures 
and small intensity these smallest grains do not contribute to 
the SED at longer wavelengths. Excluding these smallest grains, 
show indeed no variation in the SED at the maximum and longer 
wavelengths (red line). The SED of the AMM grains decrease s 
significantly at 70 /rm as observed by e.g. lStennik et al.l ll2003h . 

When adding the ice mantle the far-IR emission increases 
and the emission peak shifts to longer wavelengths. Since an ice 
band occurs at the wavelength, where the maximum is reached, 
it is not clear if there is a real shift due to a temperature decrease 
or if the ice band is dominant. The calculated temperatures of 
AMM and AMMI grains show a small decrease of less than 
1.5 K for the larger aggregates. The temperature of the small¬ 
est purely carbonaceous AMM grain decreases to 14 K when 
adding the ice mantle. Excluding these smallest AMMI grains, 
we find no variations in the SED (blue line compared to cyan 
line). 

Eitting the calculated SEDs with a modified blackbody, we 
can derive the dust colour temperature, the spectral index and 
opacity at 250 pm. In Eig.|6]and Table [1] these values are sum¬ 
marised for the SEDs of the different grain populations, inte¬ 
grated over the Herschel and Planck bands and over the 100 
to 500 pm and 100 to 870 pm wavelength ranges. The sepa¬ 
rate BGs have colour tem peratures of around 20 K, i n agree - 
ment with the findings of I Planck Collaboration et akl (l2014ah . 
while the colour temperatures of CMM grains decrease to about 
18.7 K. The SED consists of three different components, each 
having a different size distribution. The colour temperature is 
the value to fit the SED with a modified backbody and not 


the temperature of the grains. Compared to CM grains CMM 
grains with an amorphous carbon core are hotter while CMM 
grains with an amorphous silicate core have similar tempera¬ 
tures. The SED therefore does not show a shift of the peak 
position, but because of the additional change in the spectral 
slope, the best modified blackbody fit gives a lower tempera¬ 
ture. Eor the aggregates, the colour temperature is about 2 to 3 
K lower than the CM grains. Compared to the calculated tem¬ 
peratures for grains larger than 100 nm, the colour temperatures 
are on average around 1 K higher. The spectral index increases 
from around 1.5 for separate grains to about 1.8 for CMM and 
AMM grains and up to 2 for A MMI grains . The o pacity at 250 
pm for the original model of Llones et ^ (1201 3h is found to 
be 0.5 X 10^^ cm^/H, which is esse ntially the same as the ob¬ 
served value of 0.49 x 10^^ cm^/H (iPlanck Collaboration et al.l 
l2014allbl) . We see an opacity increase of around 2.3 for CMM 
grains and 2.9 for AMM grains, which is in good agreement with 
the observational results discussed in Sect. 1. Eor AMMI grains, 
the increase is up to a factor of 7, which is due to the increase in 
mass when accreting the ice mantle. 

4.3. Extinction 

Using DustEM, we also derive the extinction of the CMM grains 
and the aggregates and compare them to observations. In Eig.|2l 
we show the derived extinction compared to observations nor¬ 
malised at 10 pm (left) and at the K-band (right). The normal - 
isation at 10 pm, as suggested by Ivan Breemen et all (1201 lit , 
shows the variations at different wavelengths in the extinction 
curves between AMM and AMMI. The ice features are visible 
and an increase between 6 and 8 pm and between the two sil¬ 
icate features occurs. The usual normalisation to the K-band is 
problematic, since we find differences in the dust extinction in 
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this wavelength range when coagulating grains and/or accretion 
mantlefl However, in order to compare the modelled results 
to observations the data have to be normalised. T he compari¬ 
son to obse rvations show that th e original model bv Llones et alJ 
(1201 3h and iKohler et al.l (|2014|) agrees well with the observa¬ 
tional data for the diffuse ISM from iMathisI (Il990l) . The silicate 
featu res of the aggregates are in agreement with observations 
from [Rieke & Lebofsk\^ (Il985h but vary at waveleng t hs sho rter 
than 5 jim. The silicate features observed bv lMcClura (l2009l) are 
in general much less pronounced than for the model of the dif¬ 
fuse ISM and the aggregates. However, the differences between 
our model and the observational data might be due to the dif¬ 
ferent methods and assumptions used in the normalisation of the 
observational data. This topic will be discussed in detail in a fol¬ 
lowing paper. For our model we derive Ry values of 3.5 + 0.2 for 
CM grains, 3.9 + 0.2 for CMM grains, and 4.9 + 0.2 for AMM 
and AMMI grains. The relatively large uncertainties are due to 
the appearance of the ripple and interference structures in the 
extinction curves, which result from the rather narrow size dis¬ 
tributions. 

In Fig.[8]the extinction normalised to 1 jum“' in the visible- 
UV is shown. For evolved grains, we find strong differences in 
the UV bump and FUV extinction, which are characteristic of 
the diffuse ISM and where they appear to arise from small car¬ 
bonaceous grains (VSGs). We assume that the VSGs are com¬ 
pletely coagulated onto the surface of the BGs and this is why 
the UV bump disappears for our aggregates. Coagulation and ac¬ 
cretion leads to a strong increase in extinction compared to sep¬ 
arate grains. We find the same dec r ease a s seen in observations 
published by iFitzpatrick & Mas'sal (l2007ll . except that in those 
observations the UV bump is still slightly visible. This might be 
due to VSGs belonging to the diffuse gas surrounding denser re¬ 
gions or distributed along the line of sight and thus contributing 
to the observed extinction. 

In Fig.|9]the mass extinction coefficient k is presented for the 
original CM grains as well as for the CMM, AMM, and AMMI 
grains. For CM grains the amorphous carbons show the highest 
K values at the longest wavelengths and the lowest k values with 
the addition of the aliphatic-rich mantle. For CM grains with 
amorphous silicate cores, k increases by about a factor of 2 with 
the addition of the aliphatic-rich mantle. For AMM grains, k is 
similar to the values for CMM grains with amorphous silicate 
cores, and increases slightly with the addition of an ice mantle. 
It is therefore less than ideal to choose a single value of k to 
derive the total mass of the dust in a galaxy, since dust evolution 
leads to large differences in k. In Table |2] we present k for the 
different models. 


4.4. Radiative transfer approach 

We also combine DustEM with the CRT radiative transfer code 
as described in Sect. 3.3. The results for the derived spectral in¬ 
dex p and colour temperature are presented in Fig. [10] for the 
four different models presented in Sect. 3.3 and for the analysis 
with SEDs integrated in the Planck and Herschel bands and in 
the 100 to 500 ;um and 100 to 870 pm wavelength ranges. 

Eor the CM and CMM grains (Models CM and CMM), yS 
decreases with decreasing temperature in all cases. This cor- 

^ As pointed out bv iJones et all ( 120131) . it is likely that the extinction 
in the B and V hands, as we note for the K hand here, varies and so 
normalisation of the extinction curve at the B and V (and K) bands 
biases the interpretations. 



Fig. 8 . UV extinction of single grains and aggregates normalised 
at 1 /im“'. 



Wavelength in pm 

Fig. 9. K shown for the different grains. The grey lines indi¬ 
cate the wavelengths 250, 350 and 500 pm. 01 indicates amor¬ 
phous silicate of olivine-type, Py indicates amorphous silicate of 
pyroxene-type and Ca indicates amorphous carbon. 

relation is a pure radiative transfer effect: when the clouds 
are denser, the mixture of temperatures along the line of sight 
increases resulting in broader SEDs and thus lower p val¬ 
ues. This plot clearly shows the increase in beta, when ac¬ 
creting a second mantle of aliphatic-rich amorphous carbon. 
Assuming grain coagulation into aggregates and ice-mantle for¬ 
mation (Models AMM and AMMI), we can clearly see that 
P increases with decreasing temperature down to about 15 K. 
At even lower temperatures, p starts to decrease again. This 
means that the radiative transfer effects are strong enough to 
overcome the increase in p of the aggregates. The /I-T anti¬ 
correlation is strong for all considered bands but is more dis¬ 
tinct when analysing the data in the Planck bands, i.e. over a 
wider wavelength range. The colour temperatures decrease from 
around 19.5 to 15 K and p increases from 1.45 to 1.8. These 
values are reached in re gions with Ay 16. The results m atch 
the values measured by [ Planck Cc tllaboration et alJ (1201 lah and 
iPlanck Collaboration et al.R2014ah . 

5. Conclusions 

We investigate how evolutionary processes, i.e. accretion and co¬ 
agulation, change the optical properties of grains in the transition 
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Table 2. The mass extinction coefficient k in crc? jg at A - 250, 350 and 500 jum for the four models presented in the paper. The 
mass-weighted k is given in the column labelled Tot; it represents the k for extinction and not emission, since the dust temperatures 
for each component are not taken into account. 


A 

/im 

Ca 

CM 

01 

Py 

Tot 

Ca 

CMM 

01 Py 

Tot 

AMM 

AMMI 

250 

17.4 

5.4 

5.3 

8.7 

0.7 

10.2 

10.2 

6.4 

10.4 

12.5 

350 

11.7 

2.8 

2.6 

5.3 

0.5 

5.1 

5.1 

3.2 

5.3 

6.2 

500 

7.8 

1.3 

1.3 

3.4 

0.3 

2.5 

2.5 

1.5 

2.6 

3.1 
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Fig. 10. yS-T correlation for cloud cylinders derived with radiative transfer calculations for the four models described in Section 3.3. 


from the diffuse ISM to denser regions. For dust in the di ffuse 
ISM we assume the dust properties of the Llones et alJ (l20I3h and 
iKohler et al.l (|2014|) model. We also assume that in denser re¬ 
gions these grains undergo evolutionary processes: the accretion 
of material from the gas phase forms an additional aliphatic-rich 
mantle of amorphous carbon and the grains then coagulate into 
aggregates. In even denser regions ice mantles form on the sur¬ 
face of the aggregates. With detailed model calculations we de¬ 
rive the optical properties of these evolved grains and use the 
DustEM tool and the CRT radiative transfer code to derive the 
SED and extinction. With the presented dust-evolution model 
we are able, for the first time, to self-consistently explain the 
four observed changes in the SED in going from the diffuse ISM 
to denser regions with Ay of up to 16, namely the changes in 
temperature, spectral index (and their anti-correlation), opacity 
(emissivity) and the mid-infrared emission. 

The y6-T anti-correlations that we find in the Herschel bands 
differ from the observational results. The PACS and SPIRE chan¬ 
nels have a high level of noise which mig ht significant l y influ - 
ence the measured relation, as shown by lYsard et^ te012l) . 
This is not the case for the Planck HEI-i-IRAS bands, where 
we find a good correspondence between the model and obser¬ 
vations. We note that the changes that we observe in fi are ac¬ 
tually driven by gas density-dependent dust evolution processes 
(accretion and coagulations). Thus, the jS-T effect is actually a 
jS-hh effect according to the predictions of our model. 
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